Perspective on Afterglows: Numerically Computed Views, 
Lightcurves and the Analysis of Homogeneous and Structured 

Jets with Lateral Expansion 
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ABSTRACT 

Herein I present numerical calculations of lightcurves of homogeneous and structured after- 
glows with various lateral expansion rates as seen from any vantage point. Such calculations 
allow for direct simulation of observable quantities for complex afterglows with arbitrary energy 
distributions and lateral expansion paradigms. A simple, causal model is suggested for lateral 
expansion of the jet as it evolves; namely that the lateral expansion kinetic energy derives from 
the forward kinetic energy. As such the homogeneous jet model shows that lateral expansion is 
important at all times in the afterglow evolution and that analytical scaling laws do a poor job 
at describing the afterglow decay before and after the break. In particular, I find that lateral 
expansion does not cause a break in the lightcurve as had been predicted. A primary purpose of 
this paper is to study structured afterglows, which do a good job of reproducing global relation- 
ships and correlations in the data and thus suggest the possibility of a universal afterglow model. 
Simulations of structured jets show a general trend in which jet breaks become more pronounced 
with increasing viewing angle with respect to the jet axis. In fact, under certain conditions a 
bump can occur in the lightcurve at the jet break time. I derive scaling relations for this bump 
and suggest that it may be a source of some bumps in observed lightcurves such as that of GRB 
000301 C. A couple of lateral expansion models are tested over a range of efficiencies and viewing 
angles and it is found that lateral expansion can, in some cases, substantially sharpen the jet 
break. I show fiux surface contour maps and simulated images of the afterglows which give insight 
into how they evolve and determine their lightcurves. 

Subject headings: gamma rays: bursts — gamma rays: theory 



1. Introduction 

It is currently widely believed that gamma-ray 
bursts (GRBs) derive from narrow (half-opening 
angle Oq^grb ^ few degrees) jets of relativistic 
ejecta pointing toward the observer. One of the 
basic motivations for this has been to relieve the 
energy crisis in GRBs by reducing the necessary 
total energy from the inferred isotropic equiva- 
lent energy (up to several 10^^ ergs) by a factor 
- 10""^ - 10-^ (and thus boosting the 
total event rate by the reciprocal, to 
include unseen jets directed away from the ob- 
server). If the GRB emission is collimated, it is 
plausible that the afterglow shock is also colli- 



mated into a cone with opening angle 6>o, generally 
thought to be larger than Oq^grb- 

A rich area of inquiry is then to predict and look 
for observable consequences of a narrow jet both in 
the GRB phase (e.g. Salmonson 2000, 2001) and in 
the afterglow phase (e.g. Rhoads 1999; Sari et al. 
1999). The most fundamental consequences derive 
from the deceleration of the afterglow as it prop- 
agates into the interstellar medium (ISM). The 
relativistic motion of the emitting shock causes 
the radiation to be beamed into an angle l/F in 
the observer frame. At early times, this relativis- 
tic beaming angle is smaller than the physical jet 
opening angle, l/F < ^o, thus the emission ap- 
pears to be isotropic. Eventually, as the jet de- 
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celerates, it transitions to l/F > Oq and the fi- 
nite, non-isotropic extent of the jet becomes ap- 
parent to both the observer and to the jet itself. 
The observer sees a deficit of fiux compared to 
that expected from an isotropic emitter. The jet, 
now being entirely causally connected, begins to 
expand sideways. These two effects combine to 
cause a break in the light curve (Rhoads 1999; Sari 
et al. 1999). Roughly twenty such jet-breaks have 
been observed (see Frail et al. 2001, and references 
therein) and these constitute the most direct evi- 
dence for beaming to date. 

Also, since the later, slower emission is beamed 
into a wider angle than the earlier emission, we 
expect to see off-axis optical and radio afterglows 
where no gamma-ray or X-ray burst was seen 
(Rhoads 1997; Dalai et al. 2002). To date these 
so-called "orphan afterglows" have yet to be pos- 
itively observed. 

Thus far the study of this afterglow jet-break 
has assumed an afterglow pointed directly at the 
observer with constant, homogeneous, unstruc- 
tured energy and mass density oc H{0o — 0) across 
the jet surface, with a hard edge at ^ = where 
Hi) IS the Heaviside step function. As such, from 
scaling laws, tj oc Oq^^ ^ the observation of a jet- 
break time, tj, gives direct information about the 
opening angle of the jet. While this model is rel- 
atively simple to calculate and is amenable to an- 
alytical calculation, it is not necessarily the eas- 
iest jet morphology for nature to produce and 
thus may not accurately represent a physical jet. 
Firstly, one certainly expects the observers view- 
ing angle, 0^, of the jet to vary. Furthermore, a 
homogeneous jet with a hard edge will be unsta- 
ble to expansion and rarefaction and is thus un- 
likely to propagate intact and is unlikely to have 
been formed in the first place. Recent numerical 
simulations by Zhang et al. (2002) of relativistic 
jets emerging from stars, within the context of the 
'collapsar' model, show substantial structure, with 
the most energetic material along the jet axis and 
decreasing with larger angles from the axis. 

The first semi-analytical calculations of an af- 
terglow jet with structure, i.e. with decreasing en- 
ergy density and/or Lorentz factor as a function 
of angle for the jet axis, was done by Rossi et al. 
(2002) while and analytical treatment was done 
by Zhang & Meszaros (2002). Not long before, a 
quantitative discussion of a structured jet model 



was put forward by Salmonson & Galama (2002). 
It was found by Rossi et al. (2002) and Zhang & 
Meszaros (2002) that a universal structured jet, 
viewed at different angles, will yield a range of af- 
terglows with jet-break times relating to viewing 
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angle tj oc 6>v . By letting the energy per solid 
angle of the jet decrease with angle from the jet 
axis, ^, as e oc ^ they were able to effectively re- 
produce the observed relation E (x tj^ (Frail et al. 

2001) . 

In this paper I present calculations that fur- 
ther refine the work by Rossi et al. (2002) and 
Zhang & Meszaros (2002). By discretizing the sur- 
face of the afterglow, arbitrary functions of energy 
density, Lorentz factor and lateral Lorentz factor 
(defined in the fluid frame) can be simulated and 
lightcurves produced for arbitrary viewing angles 
0. In so doing, I corroborate some of the key 
results of (Rossi et al. 2002; Zhang & Meszaros 

2002) , i.e. tj (X and Etot oc t~^ if e oc Q-'^ . 
Also, these more detailed calculations allow for 
a quantitative discussion of how the lightcurve 
breaks; in particular I find that a flattening, and 
even a bump in the lighcurve is possible just prior 
to the jet-break time. In addition I study lateral 
expansion and find that some of the general scal- 
ing laws that are widely used in afterglow work 
are incorrect. 

2. Numerical Jet Calculations 

The calculation presented here begins by dis- 
cretizing the surface of the afterglow mapped with 
polar coordinates, (6>,^), into small elements of 
solid angle dVt = cos sin OdOdcj) where herein we 
assume the afterglow is spherical, i.e. zero inclina- 
tion -0 = 0, and where ^ = corresponds to the 
jet axis. Each surface element plows into the ISM, 
sweeping up mass according to its cross-section 
dA = R^dQ.^ decelerating, shocking and radiating. 
The physics of this calculation can be broken up 
into two parts i) the dynamical evolution of each 
surface element of the afterglow, as dictated by 
conservation of energy and momentum and ii) the 
radiative mechanism, which I take to be the stan- 
dard synchrotron shock model (Meszaros & Rees 
1997). 

To calculate the evolution of the afterglow 
shock one needs only to invoke conservation of 
energy and momentum along with an assump- 
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tion of radiation losses. Herein I assume radiative 
losses are dynamically insignificant i.e. the evolu- 
tion is adiabatic. Define the intitial bulk Lorentz 
factor To = £o/^oc^ of a surface element of the 
afterglow, where £o and Mq are the initial en- 
ergy and rest mass per solid angle. Thus following 
Paczynski & Rhoads (1993) and Rhoads (1999), 
the energy and radial momentum of a surface ele- 
ment of the afterglow is 



To + / = (1 + mr 

To(3o = (1 + f)^r(3 



(1) 



where ^ = (E' + M)/M is the internal energy of 
the expanding shell. ^ 

The mass fraction, /, of accumulated interstel- 
lar mass density, pisM = nrrip^ where rup is the 



^Eqns (1) are the correct equations for the evolution of jet. 
It has been noted by several authors (e.g. Huang et al. 
2000b) that eqns (1) yield a velocity, v oc ^ in the 
non-relativistic limit that is not consistent with the Sedov- 
Taylor (S-T) blastwave solution, v oc i?-3/2, (e.g. Shu 
1992). Thus it was assumed that these equations are in- 
correct for describing afterglow shock dynamics for the en- 
tire evolution. There have been alternative formulations 
of the afterglow shock dynamics. The reason eqns. (1) do 
not reproduce the S-T solution is because the S-T solu- 
tion (and the Blandford-McKee (B-M) solution (Blandford 
& McKee 1976)) is a "blastwave" and as such the fireball 
does work as it expands and thus radial momentum is not 
conserved. The S-T solution depends on the assumption 
of spherical symmetry (the lateral expansion of a fiuid el- 
ement is opposed by expansion of adjacent elements, thus 
forcing the fluid radially) and thus does not apply to a jet. 
The jet, having nothing to push on, expands laterally and, 
thus lacking a radial force, must conserve radial momen- 
tum. An afterglow jet is more akin to a bullet flred from a 
gun than to a spherically expanding blastwave; the jet and 
the bullet are isolated, momentum conserving bodies de- 
celerating due to their interaction with their surroundings, 
while the blastwave continues to do work as it expands, 
by virtue of its sphericity, and thus is not momentum con- 
serving. Thus I argue that eqns. (1) describe the entire 
evolution of a jet more accurately than do the S-T or B- 
M solutions in their respective regimes (so long as the jet 
hasn't expanded so far as to become spherically symmetric, 
which I do not find to be the case over the time interval 
discussed herein). Note that the degeneracy of the energy 
and momentum equations for F ^ 1 make the behavior 
of eqns. (1) indistinguishable from the B-M solution in the 
relativistic regime. 



mass of the proton, is 
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where the number of electrons swept up into the 
shock element is 



/' 

^0 



(3) 



and the the energy per solid angle is £52 = 
£/(10^^/47r ergs), and the element solid opening 
angle is 

An = sin Act) AO (4) 

where the position of a surface element, 6>, will 
evolve as the shock laterally expands due to inter- 
nal pressure (Section 3.1). The velocity magnitude 
is P = \v\/c = — 1/r^ and proper time in the 
fluid frame is 



dr 



Eqns. (1) can be solved for the Lorentz factor 

v/i + 2ro/ + p ' 



(5) 



(6) 



So by specifying £0 and Fq and a perscription for 
lateral expansion v± (Section 3.1), the entire evo- 
lution of the afterglow as a function of R is deter- 
mined. 

In order to calculate the observed flux, F^y, at a 
given frequency, z/, note that the intensity trans- 
forms as / = rS^ where the Doppler factor is 



6=[T{l-p-h)]-' 



(7) 



where (3 = v/c and h is the unit vector point- 
ing toward the observer. Following Rossi et al. 
(2002), here I focus on the power-law branch of 
the spectrum between the peak, z/^, and cooling, 
z/c5 frequencies. As such, the proper intensity at 
the proper peak frequency, z/^, is 



ji _ ^ ^ 



AnR'^An 



(8) 
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where the shock element has surface area, R^AQ, 
and the proper power per electron radiated at, z/^, 
is (Wijers & Galama 1999) 



5.4 X 10-24, _^ 



1/2 1/2 
^B,- 



(9) 



ergs s-i Hz-i electron-^ . (10) 



where the factor 0^ has been included here to 
make this equation valid in the non-relativistic 
limit (Rybicki & Lightman 1975). The observed 
intensity at a frequency, z/, is 



(11) 



and integrating over the population of radiating 
electrons, the flux will go like 



/ 



where is the luminosity distance. The surface 
area of an afterglow element as seen by the ob- 
server, AA, is calculated self consistently in the 
code by projecting the elements onto a surface 
perpendicular to the observer line of sight. The 
minimum electron frequency is 



W = 3.5 X 10^ 



0.64 



<-i^o4('_2 Hz (13) 



and ly' = ly/S where we take u = 4.4 x 10^^ Hz for 
the R-band. So 



7.96x10" 



0.64 



2 1/2 T^Sr. 



(14) 

and thus eqn. (11) scales like oc T^^+M^+^i? z/"^ 
oc r^^+^^^+^t zy", where t ex R/T/5, thus demon- 
strating the explicit scalings and consistency with 
Rossi et al. (2002). For a =1/2 the observed flux 
is 



=6.4 X 10 



-57 



0.59; V0.64 



1/2 



3/4 

^6,-1^0^5 _2 



2 Djil) 
l + zDl{z) 



(15) 



where a cosmology, {Qij^A) = (0.3,0.7), was used 
with Ho = 65 km s"^ Mpc~^ to give Dl{1) = 



2.2 X 10^^ cm. It is important to note that by 
explicitly evolving F (eqn. 6) and Fj^ (eqn. 15) in 
terms of the number of swept up electrons, TVg oc 
volume (eqn. 3), this formulation consistently ac- 
comodates sideways expansion of the jet (Section 
3.1). Finally, the flux lightcurve as a function of 
observer time, tobs, is calculated by 



Jo 



{1 + z) (l-/3-n) 



dr 



(16) 



To better compare physical timescales, the red- 
shift dependence is removed from the times plot- 
ted in this paper; t = tobs/ (1 + ^)- Thus we have 
caste the calculation of the afterglow lightcurve 
into the state variables of the problem: i^, ^7, F, 
S. 

A calculation proceeds as follows. An initial 
afterglow is specified by ^{O^cj)) and Fo(6>,^) and 
a lateral expansion perscription and is allowed to 
plow into the ISM by incrementing the radius by 
AR/R ~ a few percent. The intensity, eqn. (11) 
and observer time, eqn. (16) are saved at each sur- 
face element. Thus a lattice of (Ij^^t) pairs are 
evaluated in {0,(j),R) space. Intensity is then in- 
terpolated on dataslices of constant observer time 
t. Finally, the total flux, F^y, at each observation 
is derived from eqn. (15) where surface areas, dA, 
are calculated from a projection of the positions 
of the observed intensities onto the observer plane 
of view. 

3. Review of Homogeneous Jet Model 

It is worthwhile to begin this discussion with a 
brief review of the homogeneous jet model, which 
is amenable to analytical calculations and thus al- 
lows for comparison and validation of numerical 
results with known results. The apparent sur- 
face area of the afterglow goes like dA ^ {R Oa)'^ 
where 6a is the angular size of the effective view- 
able aperture onto the afterglow surface 



Oa 



'1/F (9^ + l/F<6>o 

(relativistic beaming dominated) 
0^ + 1/F > ^0 

(physical jet extent dominated) . 

(17) 
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Also, one can divide the Doppler factor (eqn. 7) 
into asymptotic limits 



^i/{T0l) 0, »i/r + 0o 



Using eqn. (15) the flux at a given frequency is 

and from eqn. (16) the observer time goes hke 

..A^ (20) 

Before the afterglow shock has reached its decel- 
eration radius, Rd, it coasts freely, F Fq, and 
so ^ ~ const., thus Fy ^ ^ ^ where t ^ R. 
After the shock passes the deceleration radius, i.e. 
R> Rd^ the radius and Lorentz factor are related 
by oc F~^/^ (eqns. 2,6). In this regime there 
exist asymptotic power-law slopes for F{t) only in 
the simple cases where one of the three angular 
scales, OQ^Oy^l/V^ dominates over the other two. 
These three cases are summarized in Table 1 and 
can be seen in Fig. 1. 

One characteristic timescale of the afterglow 
model is the deceleration time 

*d = ^ « 0.2 {E^2/nY'^T-%{l + Tlel) s 

(21) 

where I define Rd = (?>Eq/ AttT^pismc^Y^^ as the 
radius at which / = I/Fq (eqn. 2). For jets viewed 
well off axis, 6^J ^ 1/F, ^o, the flux (eqn. 19) varies 
like Fd oc (5^/^. Since t(xl/6, then Fd oc t-^/^. 



Table 1 : The three asymptotic limits where one of 
the three angular scales, 0{)^Oy^l/T^ dominates the 
other two. Asymptotic expressions for Oa and 5 
are given by eqns. 17 and 18 respectively. In the 
last column, the expression in parenthesis is for 
the spectral slope = 1/2 (eqn. 11). 



limit Oa ^ Flux {a = 1/2) 



(9o > , 1/F 1/F 2F (t-3/4) 

1/F > 6,, Oo Oo 2F t-3(2a+l)/4 (^-3/2) 

0, » Op, 1/F Op 2/{rol) (t^/^) 



Another key timescale is the jet-break time, tj, 
which occurs when the shock has decelerated to a 
Lorentz factor F 1/Op 

tj = iMf^hd ^ 8.5 (E52/n)V3^8/3 ^.^ 

(22) 

where ^o,i = ^o/l°- The factor 5/4 derives from 
eqn. (16) and the fact that the jet-break time oc- 
curs when 1/F ^ Op is observed at the edge of 
the jet rather than at the center. This is a factor 
of 5/2 greater than estimates derived at the jet 
center (e.g. Sari et al. 1999). For 0^ > Op the jet 
break at time tj is largely washed out, thus a more 
pertinent timescale in this regime is that at which 
the flux is a maximum. This occurs roughly when 
(9v ^ 1/r > 6>o, where S ^ 2F/(1 + O^r'^) - F so 
i^max - F^ - where t - T'^/^ 

The final phase of the afterglow is when the 
shock motion becomes non-relativistic, (3 ^ 1. 
This regime is beyond the purpose of this pa- 
per, however it is necessary to analytically de- 
scribe the non-relativistic lightcurve behavior of 
the present model so to understand the asymp- 
totic behavior of the simulations. For f3 <^ 1 
eqns. (1) yield P oc R~^, eqn. (15) goes like 
Fjy (X R^jS^. and eqn. (16) becomes t ex R/ j3. Thus 
the non-relativistic lightcurve of the present model 
is Fjy^non-rei. OC Noticc that this is indepen- 

dent of the spectral slope. This non-relativistic 
model is likely incomplete; for instance Dai & Gou 
(2001) suggest modifying the shock amplification 
of the shock-frame magnetic field from B' (x j 
to B' oc y^7(7 — 1). Such modifications can pro- 
duce substantially different behavior of the after- 
glow lightcurve, including a break at the relativis- 
tic/nonrelativistic transition (Huang et al. 2000a). 
Since our understanding of such field generation is 
still incomplete (for example Rossi & Rees 2002), 
here I choose not to make any such modifications. 
Let it suffice to understand that all post-break 
lightcurves modeled in this paper will asympotot- 
ically approach Fi, oc as they approach the 

non-relativistic regime (for example, see late times 
of the bottom plot of figs. 4, 5, 6, 7). 

In order to validate the results presented here, 
I compared with existing solutions. Comparing 
these homogeneous afterglow solutions with those 
of Granot et al. (1999), I find very good agree- 
ment when one accounts for three parameter dif- 
ferences. As mentioned above, eqn. (11) is valid 
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for observed frequencies above the sychrotron fre- 
quency, v > This simphfication is necessary to 
remove unwanted spectral breaks when attempt- 
ing to study the cause and quality of dynamical jet 
breaks. This limit corresponds to the case when 
parameter ^ ^ 1 in Granot et al. (1999). Sec- 
ondly, Granot et al. (1999) take the velocity of the 
emitting electrons behind the shock to be F/V^, 
where I have assumed it to be F. I have tested 
both choices (the former can be implemented by 
effectively replacing F F/>/2 in eqn. (15)) and 
find a quantitative difference in the flux surface, 
but the difference is minimal for the integrated 
flux; so the shape of the lightcurves is not signifi- 
cantly affected by this choice, as can be seen from 
the results presented in this section. The last dif- 
ference is the choice of spectral index; herein the 
value a = 1/2 is employed while Granot et al. 
(1999) use a = 3/4. Adopting the values used by 
Granot et al. (1999), I find very good agreement 
with their results. 



SGC min tj 


hr 






- ' ' 

r// 


1 


1 


1 
-j 
-j 






t-3/2 
















/ ^ 




\ 9^=2" 1 A 




t9/2 




n / 








r 0^=4° / 


IK 






= Q =8° // 


1 : 








t (days) 



3.1. Lateral Expansion 

Thus far the discussion has concerned jets with 
no lateral (sideways) expansion, = 0. It has 
been argued that when Cgt' ^ ^oct, i.e. when the 
Lorentz factor has decreased such that 1/F > 6>o, 
where Cg ~ c is the sound speed, the shock will be- 
gin expanding sideways (Rhoads 1999; Sari et al. 
1999). This increasing angular size of the shock, 
sweeping up an ever larger region of ISM, is sup- 
posed to impede the forward propagation and the 
shock begins to decelerate exponentially with ra- 
dius. It is argued that this shift from forward 
expansion to lateral expansion will cause an ex- 
ponential decay in Lorentz factor as a function 
of radius, V{K) ex exp(— i?) and thus a break in 
the afterglow light curve. Thus one can derive 
the predicted post-break lightcurve slope by as- 
suming forward expansion essentially ceases, R ~ 
constant. Also, if the jet is narrow, then the jet 
opening angle will expand like ^ oc cx 1/F (e.g. 
Rhoads 1999) and since this physical extent of the 
jet is growing, one has ~ 1/F (eqn. 17). Thus 
from eqn. (19) the flux goes like Fy ~ p4Q;+2 ^j^g^.^ 
J ~ F for = 0°- Also the observer time varies 
like t - F-2 (eqn. 20). Thus - ^-^^a+i) _ ^-2 
for OL = 1/2. 



Fig. 1. — Shown are calculated lightcurves for sev- 
eral viewing angles, Oy = 0°,2°,4°,8°, of an ho- 
mogeneous afterglow with opening angle Oq = 1° 
which plows into an ISM with number density n = 
1 cm~^. This afterglow is given an isotropic equiv- 
alent energy of 10^^ ergs and an initial Lorentz 
factor Fo = 1000. These curves can be compared 
with those of Rossi et al. (2002) and thus flux is 
in arbitrary units and the big tick marks corre- 
spond to decades. One can see that the analyti- 
cal, power-law solutions for each expansion epoch 
(Table 1) are well reproduced by the numerical 
calculation. Also, shown along the top axis, for 
= 0°, the theoretically predicted deceleration 
time, td ^ 0.2 s (eqn. 21) and jet-break time 
tj ^ 8.5 min (Eqn. 22) correspond well to values 
calculated numerically. Also shown are asymp- 
totic relationships for the flux at the deceleration 
time, Fd ex for 0^ ^ Oq. It is interesting 

to note that extrapolating the 6y = 0° curves 
for the free expansion epoch, F cx t^, and the 
post-jet-break epoch, F cx as well as the 

curve Fmax one finds a common point 

of intersection. This physically corresponds to the 
lightcurve of a jet of equivalent energy, but with 
an extremely narrow opening angle = 1/Fo- 
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Fig. 2. — Views of the flux contour of the afterglow 
for various points indicated in Fig. (1). The z-axis 
is the flux magnitude and the x-y plane is the view 
plane of the observer, who would be located in the 
distance to the upper left. Flux (z) and spacial 
(x and y) axes have been normalized to faciliate 
viewing. The grid lines demarcate the coordinate 
surface of the afterglow in (j)) at a resolution 
of AO = 1750, = 36O71OO. At a) the shock, 
which is poorly resolved at this resolution, has not 
reached the deleceration time (eqn. 21) and thus 
exhibits a dome-like morphology. After passing 
the deceleration time (b), the flux becomes con- 
cave at the center, as was first described by Granot 
et al. (1999). This concavity is because the on- axis 
material is observed to evolve more quickly and 
thus decelerate and cool faster than off-axis mate- 
rial. In general, larger angles from our viewing axis 
correspond to earlier times in afterglow evolution 
and thus higher Lorentz factors, F. Physically, 
the rim of maximum flux around the bowl corre- 
sponds to the angle l/F because material at larger 
angles having higher F and thus tighter beaming, 
l/F, cannot be seen, and material at smaller an- 
gles corresponding to later times has decelerated 
and thus has reduced flux. The rim thus moves 
outward as the afterglow evolves. Eventually this 
rim approaches the physical edge of the jet (c). 
As the rim passes the edge of the jet, the light 
curve breaks (d). In (e) we see two views of the 
afterglow at large viewing angle, 6^ = 8° at a time 
when most of the flux is just outside of the rim and 
thus is basically viewed edge-on by the observer. 
These calculations do not include the thickness of 
the shock (~ R/T'^) and as such one does not ex- 
pect to observe such a deep notch in a lightcurve 
at (e) of fig. (1), however, this suggests very high 
levels of polarization in orphan afterglows. 



7 



However, the assumption that the radial expan- 
sion in R effectively halts as a result of the lateral 
expansion of is an extreme limit characterized 
by a narrow opening angle, <^ 1/7, and a lateral 
expansion velocity, assumed constant, that even- 
tually exceeds the decelerating forward velocity. 
Several workers (Wei & Lu 2000; Moderski et al. 
2000; Panaitescu & Meszaros 1999) have shown 
that in practice the lateral expansion of the jet 
does not reproduce the scaling discussed above 
and cannot create a prominent, "sharp" break and 
instead tends to smooth and broaden the jet-break 
of eqn. (22). Certainly the approximation of a nar- 
row opening angle does not hold in general and we 
presently consider the behavior of v±. 

The lateral expansion is a result of internal 
pressure in the shocked medium which, in turn, 
is the result of the forward motion of the shock 
into the ISM. Thus I argue that a physically mo- 
tivated model for lateral expansion is to make the 
rate of lateral expansion a function of the rate of 
forward expansion. Thus, one does not expect the 
lateral expansion velocity to exceed the forward 
motion of the shock, from which it derives its en- 
ergy. There are two basic reasons for this. First, 
as the shock progresses into the ISM, in deceler- 
ates, thus the rate of injection of internal energy 
into the shock from bulk kinetic energy is mono- 
tonically decreasing (~ F). Furthermore the shock 
is sweeping up interstellar mass at an increasing 
rate (~ R^). As was pointed out by Moderski et al. 
(2000), this mass has no inherent lateral velocity 
component and thus new mass must be constantly 
accelerated. This is a significant source of "drag" 
on the decreasing internal pressure responsible for 
the expansion. Thus the assumption that lateral 
expansion is at all times constant is unphysical. 

In this section calculations are presented for 
afterglows with a simple lateral expansion model 
and it is found that lateral expansion serves pri- 
marily to smooth the break if it has any effect at all 
and, perhaps more importantly, provides very lit- 
tle if any steepening of the post-break light curve. 

Each shock element originally moves radially 
with a Lorentz factor 7. The relativistic jump con- 
ditions imply the shocked gas is heated by ~ 7 — 1 
and is thus driven to expand laterally via its in- 
ternal pressure. The actual lateral expansion that 
results is a complex balance between internal pres- 
sure forces of the shocked gas and the "drag" of 
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Fig. 3. — Shown are three different lateral expan- 
sion rates for an afterglow jet of initial opening an- 
gle 6>o = 1°, initial Lorentz factor 70 = 1000, and 
isotropic equivalent energy of 10^^ ergs. These are 
Rk = (no lateral expansion), R^ = 0.01, which 
corresponds to an 7^ = 3.7, and Rk = 0.1, or 
7^ = 10.5. The top figure shows the light curves 
as seen on axis, 0^ = 0°, the middle plot shows 
lightcurves viewed at 6y = 2°. Lateral expansion 
Doppler boosts the flux into larger angles, 6>, from 
the jet axis and thus on-axis observers see a re- 
duction in flux (top plot) while off-axis observers 
see a surplus of flux (middle plot). The bottom 
plot show the evolution of the jet opening angle, 
^0, with observer time. It is important to no- 
tice that none of these sideways expansion rates 
shows marked steepening of the lightcurve after 
the jet break time, tj, as has been predicted ana- 
lytically (see Sections 3.1 & 3.2). This is because 
the growth of ^0 with time (bottom plot), while 
rapidly accerating throughout the evolution, does 
not become exponential, as analytic arguments 
would indicate. 
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Fig. 4. — R-band lightcurves (top), polarization 
curves (middle) and power-law index (bottom) for 
a homogeneous jet with opening angle = 5° 
and £ = 10^^ergs/47r, To = 100. The polarization 
curves closely match those of Ghisellini & Lazzati 
(1999, Fig. 4). The first peak of the double-peaked 
polarization curves indicates the "near" edge of 
the jet coming into view, l/T > Oq — 6^, and the 
second peak corresponds to the "far" edge coming 
into view, l/F > Oq -\- 6^. Therefore comparison 
of the three figures shows an initial steepening of 
the lightcurve for a given viewing angle, ^v, cor- 
responding to the first polarization peak, and ad- 
ditional steepening occurring at the second peak. 
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Fig. 5. — The same jet as in Fig. 4, but with side- 
ways expansion (eqn. 23) Rk = 0.01, 7^ = 1.6. 
The lightcurves are qualitatively similar to those 
in Fig. 4, with the exception that the lightcurves 
break later, at around one day, due to the growth 
of ^0 prior to that time. The double break is still 
apparent for larger viewing angles, 6>^, due to the 
"near" edge of the jet becoming visible before the 
"far" edge does. 
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Fig. 6. — The same jet as Figs. 4 & 5, but with 
a sideways expansion Rk = 0.1, 7^ = 3.7. As 
in Fig. 5, the jet break happens at progressively 
later times for larger Rk due to the faster rate of 
growth of the jet opening angle, Oq. The double 
break is no longer readily apparent at larger view- 
ing angles, 6^. Thus the sideways expansion has 
"sharpened" the jet edge break. 
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Fig. 7. — The same jet as Figs. 4, 5 & 6, but 
with sideways expansion Rk = 1.0, 7^ = 10. 
One can see a signicant pre-jet-break steepen- 
ing of the lightcurve circa points a & b as per 
Section 3.2 and due to Doppler distortions. At 
later time, point d, the lateral expansion is largely 
exhausted (Section 3.1) and so the decay slope, 
In(Flux)/ ln(t) ^ —3/4, resembles that of the case 
for no lateral expansion (Fig. 4). Flux contours for 
selected points on curve 6y/6o = 0.67 are shown 
in Fig. 8. 
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Fig. 8. — Flux contours for selected points of Fig. 7 
for 6>v/^o = 0.67. See Fig. 2 for a description of 
axes. Plots a) and b) show flux rim (e.g. Fig. 2b) 
to be heavily skewed and distorted due to the 
Doppler factor of the sideways expansion; mate- 
rial near the core of the jet is seen to be brighter. 
In b) the rim expands beyond the pole of the jet. 
Plot c) shows the rim becoming less distorted as 
the jet decelerates, but still sloped, and the edge 
of the jet is visible. Plot d) shows the rim to be 
very flat and symmetrical as the jet has deceler- 
ated substantially thus sideways expansion and re- 
sultant aberration has subsided. Grid resolutions 
are A(9 = 57150 and 5750 for a) & b) and c) & 
d) respectively, with = 360790. 



Fig. 9. — What an afterglow might look like. 
Views of the flux surfaces of the previous Figure 
8 along the z-axis. The doppler brightened core is 
on the left side of the flux surface. Flux magnitude 
is measured in a normalized color scale; black is 
zero, white is the maximum. The x and y axes are 
shown in their true scale and measured in light- 
days. By comparing with Fig. 7, one can measure 
the superluminal expansion of the jet. 
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Fig. 10. — For the homogeneous jet model, the 
evolution of the jet specific kinetic energy, F — 1, 
and opening angle, 6>o, are shown for four jet ex- 
pansion parameters, Rk^ (eqn. 23). The kinetic 
energy (top) does decrease more rapidly with ra- 
dius, for large lateral expansion parameters, 
but does not become exponential in R as argued by 
Rhoads (1999). Instead it asymptotically steepens 
to the nonrelativistic limit, F — 1 oc R~^. For this 
reason the sequence of Figs. (4, 5, 6, 7) do not ex- 
hibit a progressively more pronounced break due 
to lateral broadening of the jet. The opening an- 
gle (bottom) is seen to expand rapidly once the 
jet decelerates to roughly Oq ~ 1/F as expected, 
but because lateral expansion depends on forward 
expansion, it slows down to a maximum final open- 
ing angle 6>max = + arctan(v^^). Since these 
simulations commence at R = 1 lightday (prior 
to deceleration; F ^ constant), one can see some 
early spreading, Oq ^ J {v±/R)dt ~ \og{R). While 
this particular behavior is an artifact of starting 
the simulation at = 1 lightday, it does indicate 
that strong lateral spreading can significantly re- 
arrange the initial jet morphology before it begins 
to decelerate, thus hinting at the possibility of a 
universal jet shape independent of initial condi- 
tions. 



constantly sweeping up and laterally accelerating 
fresh ISM. This balance depends intimately on the 
hydromagnetic nature of collisionless shocks and is 
not well understood as of yet. Here I simply per- 
scribe that the lateral kinetic energy of the shock, 
in its radially comoving frame, is a constant pro- 
portion of the radial kinetic energy. Thus I define 
this ratio 

7i-l 



Rk = 



7-1 



(23) 



which is a parameter measuring the shocks effi- 
ciency at converting radial kinetic energy, 7 — 1, 
into lateral kinetic energy, 7^ ~ 1^ ^ind thus can 
vary from zero to unity. 

Let us define lateral velocity as the velocity 
of the shock in the radially co- moving frame. 
This paper assumes axisymmetry of the after- 
glow so the proper coordinate velocity, 4-velocity 
and Lorentz factor in this frame are denoted re- 
spectively; v'^{0), U'Jd) = v'^{0)j'^{d), 7! = 
(1 - Let US assume the shock ex- 

pands uniformly, i.e. v±{6) oc 6 (see Sec. 4.2 for a 
relaxation of this condition), thus noting that Rk 
goes like energy, given a ratio, Rko^ at the edge of 
the jet, = this ratio at all angles, 0, will be 
determined by 



Rk 



R. 



kO 



(24) 



It is a convenient fact that since the radial velocity, 
V, and lateral velocity, v'^^ are orthogonal in the 
shock frame, the Lorentz factor of the shock in the 
lab frame is 

r = 77i . (25) 

where 7 = (1 — i'^)"^/^ and the observer frame lat- 
eral expansion velocity is then v± = v'^j^. Thus 
using eqns. (23,24,25), a given Lorentz factor, F, 
and lateral expansion efficiency, one can de- 
termine 7 and 7^. The total Lorentz factor, F, 
is the quantity that is evolved by the dynamical 
eqn. (6). The angle made by the velocity vector 
with respect to the original radial vector of the 
shock element, at 6>o, is 



/ ^1 \ 
ao = arctan — | . 



(26) 



In other words, we measure the sideways expan- 
sion with respect to the original radial vector of 
the shock element; the shock momentum does not 
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undergo a torque as expands. The angle between 
the shock element velocity and the radial vector at 
is 

a = 6>o - 6> + ceo (27) 

and the shock position, 0^ is incremented with ra- 
dius, as 



e{R) = ^0 + j a 



dR 
~R 



(28) 



There are many conceivable elaborations of this 
model. For instance, it is likley that the efficiency 
of conversion of forward motion to lateral mo- 
tion is dependent on Lorentz factor and radius. 

This model implies a maximum asymptotic 
value for the opening angle 6 which is achieved 
when the shock has decelerated to nonrelativis- 
tic velocities, at which eqns. (23,26) imply ceo ~ 
arctan(i/R^), and since a will become zero at 
large radii, then eqn. (28) gives 



(29) 



6>max = 6>o + arctan( V Rk) . 



Thus, as seen in Figs. 3&10, one does not see ex- 
ponential growth of the jet opening angle, i.e. oo 
exp(i?), and thus the analytically predicted run- 
away expansions are physically unrealizable. 

3.2. Amending the Homogeneous Jet Scal- 
ing Laws 

An interesting result of these calculations is 
that the afterglow decay slopes are altered from 
the theoretical asymptotes both before and after 
the jet-break time. In this section I describe why 
this occurs and give amended analytical asymp- 
totic decay slopes. In Fig. 3 is shown the same 
homogeneous, narrow, = 1°, jet as in Figs. 1 
& 2, but now with a range of lateral expansions. 
A key feature is that, for 0^ = 0°, greater lat- 
eral expansion is seen to steepen the lightcurve 
prior to the jet break, while not steepening the 
lightcurve after the jet break (in fact this slope 
becomes less steep). This can also be seen in the 
sequence of Figs. 4, 5, 6, 7 for a more 'realistic' jet 
with 6q = 5°. Thus we calculate behavior oppo- 
site to that predicted by Rhoads (1999); Sari et al. 
(1999). 

To explain this, begin by ignoring Doppler fac- 
tors, (5 oc r, thus from eqn. (19) the flux varies 



hke 

oc T^^^+^^RHI 
and the time (eqn. 20) 

R 

toe— . 



(30) 



(31) 



Now parameterizing the Lorentz factor depen- 
dence as 

r oc (32) 

where A is unity for the case of no lateral expan- 
sion and will be estimated otherwise. These equa- 
tions with eqn. (17) imply 



oc 



I 1 + 1/(3^) 

i+i/(r - 



t 



pre-jet-break: 6^/^ oc l/F 

t(m~ post-jet-break: 6>a oc 6>o . 

(33) 

In order to estimate the parameter, note from 
eqn. (2) that 



dR 



and from eqns. (28,26) that 



e{R) ^00 + 



r dR 
J iR 



(34) 



(35) 



A crucial point here is that the lab-frame angle 
between the velocity and radial vectors (eqn. 26), 
QfQ oc 1/7, and not oc l/F. The reason is 
that maximal sideways expansion occurs when 
the shock frame lateral expansion energy, 7^, is 
roughly the same as the radial expansion energy, 
7. From eqn. (25), this implies that total spe- 
ciflc energy of the shock, F ~ 7^. So if an ex- 
pansion angle of ao ~ l/F is desired, this implies 
7^ ^ F and the total specific energy of the shock 
is ~ F^ which is energetically untenable. The lat- 
ter, erroneous choice of 0^0 going like l/F results in 
the standard exponential decay of F with R which 
was discussed in the first paragraph of Section 3.1. 
From eqn. (6) 

2/ 



F^ 



(36) 



so 1/7 (X l/\/r (X /^/^. Thus one can solve 
eqns. (34,35) for the asymptotic behaviors. In the 
limit of no lateral expansion, 0[R) 6>o, and so 
eqn. (34) becomes f oc R^ and thus F ex R~^^^ so 
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yi = 1 as expected. For the limiting case when lat- 
eral expansion dominates, eqn. (35) gives 0{R) 
J dR/^/R (X / f^^^dR/R and eqn. (34) becomes 
f (X R^ and thus T ex R'^ so A = 2. Therefore, 
from eqn. (33), for a = 1/2 the pre-break light- 
curve is steepened from oc for negligible 
lateral expansion, to Fi^ oc t~^^^ for maximal lat- 
eral expansion. This range nicely brackets the pre- 
break decays of simulated lightcurves over a range 
of lateral expansion rates as seen in the bottom 
panel of Figs. (4, 5, 6, 7). It is also worth not- 
ing that inspection of Figs. (7, 8) indicates that, 
not surprisingly, Doppler effects become impor- 
tant for large lateral expansions. Such effects were 
neglected in the derivation of these scaling laws. 

The post-break afterglow decays less steeply 
than analytical models predicted. This is because, 
as discussed in Section 3.1, a realistic afterglow jet 
will expand laterally only in proportion to its for- 
ward expansion, thus the lateral expansion rate 
decreases at late time (see Fig. 10) and the post- 
break decay curve will asymptotically move to- 
ward the limit of no lateral expansion; A = 1. 

To summarize, I have implemented a simple 
model for jet dynamics in which the lateral ex- 
pansion of the jet is an effect of forward expan- 
sion. I argue that this basic cause/effect relation- 
ship is a necessary component of any model for 
jet dynamics. Resulting jet simulations show that 
lateral expansion i) smooths the break in the jet 
lightcurve, ii) steepens the pre-break light curve, 
but does not significantly steepen the post-break 
slope of the light-curve, in) can significantly al- 
ter the shape of a lightcurve. Thus I suggest that 
the oft-cited scaling laws described at the begin- 
ning of this section do not accurately represent the 
evolution of a laterally expanding jet. To properly 
diagnose lightcurve slopes and breaks, simulations 
such as those shown here are necessary. A pre- 
liminary conclusion that can be drawn from the 
simulations shown in this section is that large effi- 
ciencies in conversion of forward kinetic energy to 
lateral kinetic energy, i.e. Rk ^ 0.1, are probably 
not physical, in the context of the homogeneous 
jet model, because the lightcurves they produce 
bear little resemlance to observed lightcurves. In 
this way one can begin to constrain the physics of 
afterglow shocks with observations. 



4. The Structured, Universal Jet 

The primary purpose of the numerical frame- 
work discussed thus far is to quantitatively study 
the afterglow lightcurves from structured jets, 
i.e. with a non-uniform energy or velocity dis- 
tribution. As suggested by Rossi et al. (2002) and 
Zhang & Meszaros (2002), an asymptotic decrease 
in energy per solid angle like E{0) cx 0~'^ can repro- 
duce the observed relation E oc tj^ (Frail et al. 
2001). Furthermore, allowing £ ^ constant for 
angles within a core, 6 < Oc provides a natural 
explanation for the dearth of small angle jets, less 
than 3° also reported by Frail et al. (2001). As 
such I define an ansatz energy profile 

c -|q52 

How does one choose an initial ro(^)? based 
on simulations by Zhang et al. (2002), Fq ^ 
const, and so this choice will be used for the struc- 
tured jet runs in this paper. However, it is impor- 
tant to realize, and the following simulations con- 
firm this, that initial Fq only determines the de- 
celeration time of the shock, but subsequent evo- 
lution is independent of it. To see this, noting 
the global energy has dependence E ocT'^R^ from 
eqns. (1), then for any radius greater than the de- 
celeration radius, R > Rd^ the Lorentz factor is 
determined by F oc independentof initial Fq. 

In Fig. 12 is shown a series of lightcurves from a 
universal structured jet with isotropic equivalent 
energy £ = 10^^ ergs/47r, initial Lorentz factor 
Fo(6>) = 100, core angle Oc = 3° and Oq = 30°. 
These parameters generally reflect those inferred 
from burst observations by Panaitescu & Kumar 
(2002). This figure demonstrates the characteris- 
tics of a lightcurve derived from the structured jet 
described by Rossi et al. (2002). 

Using the scaling of eqns. (19, 20) with £ oc 
7^i?^ one finds 

oc F3^-^^^+3(9i£ (38) 

and 

t oc E'/^T-^/^S-' (39) 

where the asymptotic limit of Oa is given in Table 
1 dindS ^ 2F/(1 -h O^r^) - F for 6> < 1/F. Using 
eqn. (39) and £ oc {Oc/0^)~'^ for a structured jet, 
at an early time (F ^ 1/^v) the Lorentz factor 
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scales with observed angle as F cx 



-1/4 



and so 



oc 



-(a+2) 



SO the flux at a given time 



depends on viewing angle like 



-^i^.simulta 



OC I 



OC 



-5/2 



(40) 



for a = 1/2, and q = 2. 

There can be seen a gradual flattening, and 
eventual appearance of a bump, in the light curves 
of Fig. 12 with increasing viewing angle 0^. This 
can be understood when one considers that the 
early phase of the light curve, i.e. T > 1/0^ /is dom- 
inated by emission primarily along the line of sight 
(los) with energy £ oc {Oc/0^)~'^, while the flux at 
the break is dominated by the energetic core mov- 
ing at angle 6>v with respect to the observer. Here 
I show that these two components have distinct 
laws for their breaks, and this diff'erence between 
these two laws makes each component distinct for 
large viewing angles, Oy, thus creating a bump. 

The line of sight component can be modeled as 
a homogeneous jet with opening angle 0^. As such 
one expects a break when F ~ S ^ I/O a ^ 1/^v, 

so Fiyios oc F^"£ oc ^v*^^"^^^- Furthermore tios oc 
£l/3r-8/3 oc ol So 



',los 



OC t 



-l/2(4a+q) 
los 



OC t 



los 



(41) 



where = 1/2 and q = 2. The jet core component 
can be modeled as a narrow jet seen far off- axis; 
Ov ^ Oa where Oa = 0^ One expects a maximum 
flux when F - l/6>v, thus F^,core oc F'^^+^S oc 
0-{AoL+q) ^ where £ = 1 for the core. It follows 
that tcore OC £V3r-8/3 qc 0^^ . Thus the core is 
brightest at 



F oc t-3/8(4a+g) ^-3/2 

^ u.core ^ ^core ^ ^core 



(42) 



for a = 1/2 and q = 2. This relation traces the 
jet break in the lightcurve, Fj,{t)^ over a range 
of viewing angles. The different decay laws for 
eqns. (41, 42) demonstrate why both components 
of a lightcurve from a structured jet will diverge 
for large break times. Seen another way: 



oc 



(43) 



thus the contribution to the lightcurve from the 
core happens progressively later than that of the 
line-of-sight material. 



4.0.1. The "Bump'' 

As discussed above, under certain conditions 
a "bump" can appear in the lightcurve at the 
jet-break time. This typically happens at large 
viewing angles and small jet cores. For example. 
Fig. 11 shows three lightcurves from identical af- 
terglows except for variation in the size of the jet 
core, Oc^ demonstrating that a narrower core pro- 
duces a more prominent bump. This bump can 
be explained by considering how the flux varies as 
the observer sees an ever-increasing (~ I/7) angu- 
lar region of a surface with varying energy density. 
It can be shown that the observed flux varies like 
F oc 7^8 where p = 2 for I/7 < ^0 (isotropic) and 
p = 4 otherwise (jet-like). Also, the energy per 
steradian of the afterglow goes like e oc iOc/Oy 
where in this case q = 2. Thus the flux will vary 
like 



F(l/7) oc ( i 



(44) 



so that early in the evolution, when I/7 <C ^, the 
flux is down by {Oc/Oy as expected, but when the 
beaming angle has expanded to reach the core, i.e. 
1/7 = ^—^c 5 the flux averaged over the surface will 
be oc {6c/0y~'^. The flux of eqn. (44) decreases 
as 7 decreases until I/7 = p/{p -\- q — 1)0. So the 
smaller g/p is, the farther out the bump will occur. 
The condition for no bump to be observed is when 
the jet core becomes visible, i.e. I/7 = ^ — ^c, the 
flux has not yet begun to increase, thus 



p^q-l 



(45) 



So for p = 2 and g = 2 we have Oc > 0/3 and so a 
bump will be visible if the jet is viewed at angles, 0, 
much in excess of 30 ^ This analysis is born out in 
the simulations and gives constraints on the range 
of allowed viewing angles, ^, with respect to the 
size of jet core, Oc, and the steepness of the jet 
decay structure, q. 

There are several examples of bumps in light- 
curves including GRB 970228, GRB 970508, GRB 
980326 (for a discussion, see Zhang & Meszaros 
2002). While in general gamma-ray burst after- 
glows do not exhibit bumps in their lightcurves 
just prior to the break, GRB 000301C exhibited 
a prominent, achromatic bump at the jet-break 
time which has been interpreted as a gravitational 
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lensing event (Garnavich et al. 2000) and alter- 
natively as continuous energy injection by a mil- 
lisecond pulsar (Zhang & Meszaros 2001). In light 
of the present calculations, it is possible that the 
bump in the lightcurve of GRB 00030 IC was due 
to a simple perspective effect onto a jet with a nar- 
row core or a steep decay curve. This explanation 
is appealing in that it does not require external 
mediators (i.e. a lensing body or a pulsar) to cre- 
ate it. If the bump in GRB 00030 IC and possibly 
those of other burst lightcurves could be positively 
attributed to perspective onto a structured after- 
glow, key information could be determined about 
the size (i.e. the jet core Oc) and shape (the decay 
structure q) of the afterglow. 

4.1. Uniform Lateral Expansion 

Having constructed a structured jet model, we 
are interested in the effects of lateral expansion on 
said model. Since the evolution of the afterglow 
shock is uncertain, I choose to study two basic 
paradigms. First, in this section I simply apply 
uniform expansion to the afterglow as was done 
for the homogeneous jet in the previous section. 
Thus expansion is governed by eqn. (24). The sec- 
ond model, discussed in the next section, employs 
nonuniform expansion, where the hotter core ex- 
pands faster than the cooler wings of the afterglow. 

Examination of the progression of increasing 
lateral expansions Figs. 12, 13 & 15 shows that 
lateral expansion suppresses the bump at the jet 
break time viewed at large Oy. However, this also 
makes the breaks at small 0^ less pronounced. 
This behavior, where the jet break is more pro- 
nounced at high viewing angles, 6>v, than at small 
ones, is quite general. Even a nonuniform ex- 
pansion paradigm discussed in the next section 
demonstrates this. 

4.2. Nonuniform Lateral Expansion 

Uniform expansion of the structured jet is likely 
oversimplified. A more accurate sideways expan- 
sion perscription should encapsulate one's hydro- 
dynamic intuition that a fluid element will be ac- 
celerated proportionally to the gradient of its in- 
ternal energy density. A simple model with this 
characteristic is 



using eqn. (37) and where Oc is not constant, but 
allowed to laterally expand with the jet by stay- 
ing assigned to a particular fluid element. Notice 
that eqn. (46) replicates the uniform expansion of 
eqn. (24) for small angles, <^ Oc, but transitions 
to a constant, rigid-body, expansion for large an- 
gles. Figs. 16 & 17 demonstrate that non-uniform 
expansion can be very effective at creating a very 
sharp break in the lightcurve. 

A key approximation of eqn. (46) is that it as- 
sumes a fixed functional form for Rk even as the 
shock surface evolves. This represents a model in 
which the proper hydrodynamic timescale, which 
increases with time, become longer than the decel- 
eration timescale, thus the early time expansion 
morphology is "frozen in" and predominantly de- 
termines the subsequent evolution. Further work 
is required to determine the hydro dynamical evo- 
lution of the afterglow. 

5. Discussion 

I have implemented a simple model for after- 
glow jet spreading in which lateral expansion de- 
pends upon energy from forward expansion. This 
model demonstrates that a homogeneous jet does 
not exhibit a lateral expansion dominated phase, 
as described by analytical arguments (Rhoads 
1999; Sari et al. 1999). In particular, dynami- 
cally 7 oo exp(— i?), (Fig. 10) and observationally 
oo t-P after the break, (Figs. 5, 6, 7). The 
only source of a break in this model is the obser- 
vation of the physical edge of the jet. As shown in 
Figs. 5, 6, 7, this break is 'sharpened' by lateral 
expansion in that it becomes less dependent on 
viewing angle. 

Also, I have studied structured jets as seen at 
various viewing angles and with various sideways 
expansions. This study confirms the key results 
of Rossi et al. (2002); Zhang & Meszaros (2002), 
tj oc Oy^ and Eiso oc tj . Furthermore, the jet 
break is seen to be sharp and can exhitibit some 
flattening and even a bump prior to the break at 
large viewing angles compared to the core size, 
Oy ^ Oc. It appears to be a general feature of 
the universal jet model that late time jet breaks 
(i.e. large ^v) are more pronounced than early time 
breaks. That is, if early time breaks are sharp, 
then late time breaks will tend to have a flatten- 
ing or a bump (e.g. Fig. 13), but if this bump 
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is quenched, perhaps by nonuniform spreading of 
the jet, the early-time breaks become weak or are 
washed out altogether (e.g. Fig. 17). In general, 
the sharpness of the structured jet break is de- 
termined by the size of the core, t oc Ol^^, and 
is thus intrinsically sharper than the break in the 
homogeneous jet model which depends upon the 
outside edge t cx 0^^^. As such the structured jet 
can explain sharper breaks than the homogeneous 
jet model. 

Future work with E. Rossi et al. will focus on 
polarization (e.g. Fig. 4) as a tool to discrimi- 
nate between the homogeneous and structured jet 
paradigms. Also, improved understanding of the 
evolution of the jet, whether it be hydrodynamic 
or otherwise, will allow for more quantitative ex- 
pansion models and more predictive simulations. 

I wish to convey my sincerest thanks to E. 
Ramirez-Ruiz, E. Rossi, D. Lazzatti and M. Rees 
for their kind hospitality at the U. of Cambridge 
and insightful, enjoyable discussions during the 
preparation of this document. I also wish to thank 
T. Galama and R. Sari for stimulating, useful dis- 
cussions regarding this work. This work was per- 
formed under the auspices of the U.S. Department 
of Energy by University of California Lawrence 
Livermore National Laboratory under contract W- 
7405-ENG-48. 
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Fig. 11. — The bump in the hghtcurves at the jet 
break time becomes more prominent with smaher 
cores. These afterglows have Eiso = 10^^ ergs, 
7o = 100, v± = 0, and are seen at 6^ = 12°. 
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Fig. 13. — Lightcurves for the same structured 
jet as Fig. 12, but with moderate uniform lateral 
expansion, R}^ = 0.1. Notice a slight suppression 
of the bump at the jet-break time for large ^v/^o 
compared to Fig. 12. Fig. 14 shows selected flux 
contours for 0^ = 12°. 



Fig. 12. — Lightcurves for a structured jet with 
isotropic equivalent energy 10^^ ergs, £52 = 10, 
using eqn. (37), and initial Lorentz factor Fq = 100 
everywhere. One sees corroboration of the basic 
thesis of Rossi et al. (2002), with some flattening 
at larger viewing angles. 



18 




Fig. 14. — Flux contours for selected points of 
Fig. 13 for 6>v = 12°. At early times, a), the flux 
surface is distorted from a symmetric bowl (see 
Fig. 1 b) by both lateral expansion, as in Fig. 8 
a), and the intrinsic structure of the jet. Near 
the jet-break time plots b) show different views 
of the flux surface. Plot b3) shows the physi- 
cal pole of the afterglow nearly coincident with 
the flux peak. This demonstrates the origin of 
the break in the lightcurve in the structured jet 
model: the coincidence of the peak flux rim and 
the actual pole of the jet. As such the edge of 
the jet does not play a role in the jet break as 
it does in the homogeneous jets. Plots c) and 
d) show this coincidence clearly. Grid resolutions 
are AO = 307300, = 3607360 for a) and 
AO = 307100, = 360760 for b), c) and d). 
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Fig. 15. — Lightcurves for the structured jet 
of Fig. 12 with large, uniform lateral expan- 
sion, Rk = I. Lateral expansion suppresses and 
smooths the bump in the lightcurve at the jet 
break time for large viewing angles, 
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Fig. 16. — Lightcurves for a structured jet with 
a nonuniform lateral expansion (eqn. 46). Other 
parameters are the same as in Fig. 12. One can see 
that even a rather small, Rk = 0.01, nonuniform 
lateral expansion can affect the lightcurves when 
compared to Rk = in Fig. 12. 
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Fig. 17. — The modest nonuniform expansion 
shown here, Rk = 0.1, compared to Figs. 16 & 
12 makes a very sharp break with no bump for 
6>v > 0.2^0, but washes out the break as seen at 
smaher viewing angles. Fig. 18 shows selected flux 
contours for 0^ = 12°. 




Fig. 18. — Selected flux surfaces for a structured 
jet with a nonuniform expansion paradigm at 0^ = 
12°. Grid resolutions are AO = 30°/100 and = 
360°/120. Prior to the break, a), the morphology 
is similar to the uniform expansion (Fig. 14). The 
hot core of the jet expans and decelerates more 
rapidly than the wings. Thus a crater appears as 
the core comes into view (b - e). This sharp flux 
deficit makes a sharp jet-break (Fig. 17). 
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